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Infrared images of the dark cloud core B59 were obtained with the Spitzer 
Space Telescope as part of the "Cores to Disks" Legacy Science project. Photom- 
etry from 3.6-70 jum indicates at least 20 candidate low-mass young stars near the 
core, more than doubling the previously known population. Out of this group, 13 
are located within ~ 0.1 pc in projection of the molecular gas peak, where a new 
embedded source is detected. Spectral energy distributions span the range from 
small excesses above photospheric levels to rising in the mid-infrared. One other 
embedded object, probably associated with the millimeter source B59-MMS1, 
with a bolometric luminosity Lboi ~ 2 L , has extended structure at 3.6 and 4.5 
jum, possibly tracing the edges of an outflow cavity. The measured extinction 
through the central part of the core is Ay > 45 mag. The B59 core is producing 
young stars with a high efficiency 

Subject headings: Stars: Pre-main-sequence - Infrared: Stars - ISM: Clouds - 
Infrared: ISM 



1. Introduction 



One component of the Spitzer Space Telescope "Cores to Disks" (c2d) Legacy Science 
program (Evans et al. 2003), is imaging of nearby dark cloud cores at infrared wavelengths 
ranging from 3.6 to 70 fim. The program takes advantage of Spitzer 's high sensitivity (Werner 
et al. 2004) to survey young stars down to ~ 0.001 L in these cores. Five large molecular 
clouds are also being studied by c2d. This will allow comparisons to be made between 
the young star populations in isolated cores and large molecular clouds. Spitzer 's high 
sensitivity in the mid-infrared makes it ideally suited to studying the intermediate stages of 
star- format ion, particularly in dense cores, where dust extinction can obscure the light at 
visible and near-infrared wavelengths. 

This paper presents highlights of the c2d observational data of the dark cloud core B59. 
Though it is a known site of star-formation, there have been relatively few published studies 
of this core to date and spectral information is apparently lacking for nearly all of the objects 
discussed here. The Spitzer data show the core to be the site of a small cluster of young 
low-mass stars at a range of evolutionary stages. 
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2. B59 



B59 (Barnard 1927; Lynds 1962) is an irregularly-shaped dark cloud roughly 30' x 
15' centered around a(J2000) = 17 h 11.4 m and 5(J2000)= -27° 26' (Herbig 2005). It sits at 
the NW end (in galactic coordinates) of the Pipe Nebula, a set of filamentary dark clouds 
located close to our line of sight to the Galactic Center and above the galactic plane (/ = 357°, 
b = 7.2° for B59). The nebula also lies close in projection to the Ophiuchus molecular cloud. 

The Pipe Nebula was mapped in the J = 1 — lines of 12 CO, 13 CO, and C 18 with the 
NANTEN telescope using a 2.7' beam (Onishi et al. 1999). They identified 14 C 18 cores. 
Molecular outflow emission was detected only towards a single core (their Core 1), within 
B59. This core has roughly twice the column density of the other cores. The outflow was 
not clearly identifiable as bipolar, but may be due to two or more different outflows. We will 
refer to this as the B59 core, and sometimes loosely B59, in this paper. 

The distance to B59 is uncertain. The radial velocity of the molecular gas is consistent 
with a distance similar to the Ophiuchus cloud (Onishi et al. 1999), which is approximately 
d = 125 ± 25 pc (de Geus et al. 1989). The distance to the Pipe Nebula was recently 
estimated by combining extinctions with parallaxes (Lombardi et al. 2006). They obtained 
d = 130^20 P c > where the uncertainties are la. This distance is adopted in this paper. 

B59 is the only confirmed star-forming region in the Pipe Nebula to date. In addition 
to the CO outflow(s), there are at least 5 Ha emission-line stars in or near the core: LkHa 
346 (NW and SE), KW 2, B59-2, and LkHa 347 (Cohen & Kuhi 1979; Herbig & Bell 1988; 
Kohoutek & Wehmeyer 2003; Herbig 2005). Two other Ha emitters, KK Oph and V359 
Oph, are farther from the core center, ~ 20' and ~ 36', or 0.8 and 1.4 pc in projection, but 
may also be associated with B59 (Herbig 2005). 

There are 3 bright IRAS sources in the Point Source Catalog within ~ 5.5' of the 
center of the core with spectral energy distributions consistent with young stellar objects, 
according to the criteria of Lee & Myers (1999): IRAS 17079-2719 (associated with LkHa 
346), IRAS17081-2721, and IRAS 17082-2724. Because of the proximity of the emission-line 
stars and IRAS sources, B59 is categorized by the c2d team as a "starred" core. A 1.3 mm 
continuum survey revealed an additional probable protostellar object, B59-MMS1 (Reipurth 
et al. 1996), not associated with any of the IRAS sources. 

There have been several studies of B59 that have focused on the multiple stars. B59- 
1 and the emission-line star, LkHa 346, are probable triple systems (Chelli et al. 1995; 
Koresko 2002). Another emission-line star, B59-2, is at least double (Reipurth & Zinnecker 
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1993; Chelli et al. 1995), and may be triple (Koresko 2002). 

3. Observations 



B59 was observed with the Spitzer Infrared Array Camera, IRAC (Fazio et al. 2004), on 
UT 2004 September 3.9 and 4.2 under Astronomical Observation Requests (AORs) 5131776 
and 5132288. The pixel size is 1.2". Observations at two different epochs allowed identifica- 
tion of asteroids. The IRAC 3.6 /im and 5.8 /im arrays observe the same field simultaneously, 
and the 4.5 /im and 8.0 /im arrays a different field. The observations at each epoch con- 
sisted of a 3 x 4 map with small overlap regions. For each map there were two dithers 
of 12-sec exposures, preceded by short high-dynamic-range 0.6-sec exposures (HDR). The 
HDR observations allowed photometry of sources that saturated at longer integration times. 
The typical effective observation time was 48 sec. The point-source flux limits for reliable 
detection were roughly 0.03, 0.03, 0.10, 0.15 mJy at 3.6, 4.5, 5.8, 8.0 /mi, respectively. 

Observations of B59 at 24 and 70 /im using the Multiband Infrared Photometer for 
Spitzer, MIPS (Rieke et al. 2004), were obtained on UT 2005 April 6.4 and 7.0 (AORs 
9409280 and 9438720). The pixel size is 2.5" at 24 /im and 10" at 70 /im. The observations 
were 3x3 photometry maps with 1 cycle of 3 sec observations at 24 /im, and 3 cycles of 3 
sec at 70 /im. The second epoch was shifted by roughly 6', to increase coverage at 70 /im. 
Typical observation times were 96 sec at 24 /tm and 100 sec at 70 /im. The point-source flux 
limits for reliable detection were roughly 0.5 mJy at 24 /im and 100 mJy at 70 /im. 

The data presented here are from the Spitzer Science Center's (SSC) pipeline version 
Sll.4.0. The data were further processed by the c2d team to correct artifacts, where possible. 
This was followed by creation of new mosaics, source photometry up to 24 /im, bandmerg- 
ing for the detected wavelengths, and preliminary source typing with extinction estimates. 
Counterparts in the 2MASS catalog were assigned if found. The procedures are described 
in Evans et al. (2005). Photometry at 70 /im was done separately with the SSC's Mopex 
package, release version 030106 (Makavoz & Marleau 2005). 

One exception to these procedures is that three sources in the B59 field saturated at 24 
/im (1, 7, and 11 in Tables 1 and 2) and one source saturated at 70 /im (source 11). For 
these sources, photometry was done by fitting the radial profile in the wings with a PSF 
obtained from an unsaturated bright point source in the field. The flux uncertainty due to 
this fitting is estimated to be 20 %. 
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The Spitzer photometry is summarized in Table 2. The uncertainties given do not 
include absolute calibration uncertainties. These are currently estimated by the SSC to be 
10 % for IRAC, 10 % for MIPS 24 /im, and 20 % for MIPS 70 /im. 

The results in this paper are limited to the region where all IRAC bands and the MIPS 
24 /im areas overlap, about 170 arcmin 2 . The MIPS 70 coverage is a subset of this, about 
60 % of the area. 



4. Results 



An overview of the young star cluster in B59 is provided by the IRAC 3-color image 
and the MIPS 24 /im image in Fig. la and Fig. lb. Source 10 (Fig. lc) is near the peaks in 
the molecular gas and in the dust extinction, which are discussed below. 

Selection of candidate young stars in B59 was based on infrared excesses relative to 
photospheric colors, after correction for extinction. The excess method is needed because of 
the absence of any systematic study of the core membership. Correction for extinction is 
vital due to the high dust column densities toward the core. 

As part of its bandmerging, c2d checks whether a reddened stellar photosphere can 
fit the 2MASS and Spitzer fluxes at all wavelengths available for each source, within the 
uncertainties. When this was successful, the source was dropped from consideration as a 
candidate young star here. Note that some of these may in fact be young stars, though with 
no detected infrared excesses. 

For the remaining sources, which could be young stars, main-sequence or post-main- 
sequence stars with excesses, or galaxies, we imposed a 24 /im flux lower limit, Fj,(24 /im) 
> 3 mJy This tends to cut out galaxies as the number of galaxies brighter than this is 
~ 50 deg~ 2 (Marleau et al. 2004; Papovich et al. 2004), or ~ 2 in our field. (Two of our 
candidates have no measured 24 /im fluxes due to confusion with nearby sources, but their 
expected 24 /im fluxes would easily make this cut.) 

The galactic background counts toward B59 are less certain. The model of Wainscoat 
et al. (1992) indicates that for F^(24 /im) > 3 mJy, the number of non-photospheric galactic 
sources, primarily O-rich and C-rich asymptotic giant branch (AGB) stars, should be small, 
~ 13 deg -2 , or ~ 0.6 in our field. But models can serve only as a guide for small areas. 

Lombardi et al. (2006) identified a population of AGB stars behind the Pipe Nebula 
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from their location in the 2MASS color-color diagram. These are bright objects, peaking 
at K ~ 7, most likely in the galactic bulge. If, as suggested by Lombardi et al., they are 
similar to the OH/IR star sample of Jimenez-Esteban et al. (2005), then they would be easily 
detectable at 24 /im, F„(24 jum) > 1 Jy. But they are too few in number to significantly 
contaminate the young star candidate sample; at b = 7.2°, there are ~ 3 deg -2 , or ~ 0.14 in 
our field. 

To be conservative, we also cut out sources with excesses which appeared close to pho- 
tospheric levels. We required at least one flux to be a factor ~ 4 higher than the expected 
photosphere at that wavelength. This eliminated 4 more objects. In doing so, we may have 
dropped some young stars with weak dust excesses. The motivation is that c2d is finding that 
there are some sources with small IR excesses toward the Serpens cloud that are background 
AGB stars. 

The 20 candidate young stars are listed in Tables 1 and 2. Positions are shown in Fig. 
lc. Note that 13 of the objects lie within ~ 0.1 pc in projection of the molecular gas peak 
(see Fig. 1). We cannot completely eliminate the possibility that some of the proposed young 
stars are background objects without further data. So we refer to the set as candidate young 
stars. But some have Ha in emission (Section 2), and the majority are clustered together 
close to the center of the core. 

For the candidate young stars, we re-evaluated their extinction taking into account the 
possibility of excess emission in the near-IR. There is no systematic visible photometry of 
the cloud members, nor is there information on spectral types or effective temperatures for 
most of the objects. We rely on typical dereddened near-IR colors of young stars. 

We assumed that the young stars with all 3 2MASS fluxes in our sample have intrinsic 
JHK S colors that fall on the locus of the sample of T Tauri stars in Taurus, whose dereddened 
colors were estimated by Meyer et al. (1997), see Fig. 2. This requires that we are sampling 
similar spectral types and dust properties as in Taurus. It is difficult to test this without 
spectroscopic data on the stars. But the assumption has been used for other young clusters, 
e.g. IC 348 (Muench et al. 2003) and NGC2024 (Haisch et al. 2001), in the absence of 
better constraints. And Luhman (2004) shows consistency of extinction estimates done this 
way compared to those obtained from visible spectra for young stars in Cham I. 

We adopted the dust extinction law of Indebetouw et al. (2005), which extends from 
1.25-8.0 fim. The 1.25-2.2 /im extinction of the Pipe Nebula agrees with that of Indebetouw 
et al. within the uncertainties (Lombardi et al. 2006). Using this dust extinction law, we 
allowed the young stars to deredden to the locus of JHK S colors of the T Tauri stars, or 
to the nearest extrapolated point (Fig. 2). Uncertainties from both the data and the locus 
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envelope are propagated to the extinction estimate (Table 1). By convention, the extinction 
is given by Ay in magnitudes, whereas we do not know the extinction law in this region 
down to visible wavelengths. We only estimate near-infrared values of the extinction, e.g. 
Ax- We converted using Ay/A^ = 8.9 from the Rieke & Lebofsky (1985) extinction law, 
but note that our actual estimates are Ak- The adopted Indebetouw et al. extinction law 
appears appropriate for B59 in the JHK S bands, as the bright stars in the field follow it (Fig. 
2). 

For source 8, with only HK S data, we deredden to the center of the T Tauri locus, and 
assign a larger uncertainty. The embedded objects (10, 11) have insufficient near-IR data to 
correct for extinction by this method. 

The dereddened spectral energy distributions (SEDs) are shown in Fig. 3. The 24 jum 
extinction, A 2 4/A^ = 0.24, was taken from the Weingartner & Draine (2001) model which 
fits the Indebetouw et al. points (Ry=5.5, Case B). 

The 2-10 /im or 2-20 /im slope is commonly used to classify young stars in the IR. 
Following earlier authors (Lada 1987; Wilking et al. 2001; Greene et al. 1994) we define 
the spectral index as a = dlog(XF\)/dlog\. We define the IR spectral classes as follows: 
a > 0.3 is a Class I source, 0.3 > a > —0.3 a "Flat" type, and a < —0.3 a Class II. These 
agree with Greene et al. (1994), except that the Class III designation is restricted to sources 
with practically photospheric slopes after dereddening, a ~ — 3 (Kenyon & Hartmann 1995). 
Because of our selection criteria, we have no Class III sources in the present sample. We 
use the 2MASS K s and IRAC 8 /im fluxes as endpoints, unless noted in Table 2. Slopes 
and IR class are given in Table 2 for the observed fluxes uncorrected for extinction and with 
correction for extinction. (Source 1 after correction is on the border of Flat and Class II, 
but the complete SED indicates it is most like the Flat objects.) 

The similarity of the dereddened SEDs to known young stars in other star-forming 
regions gives us confidence that the Ay values are reasonable and that all of the objects are 
young stars. Note in particular (Fig. 3) the strong similarity of most of the Class II sources 
with a > —1.6 to the median Taurus Class II SED from the sample defined by D'Alessio et 
al. (1999). 

4.1. Extinction Map 



There are sufficient numbers of background stars detected through the B59 region to 
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provide an approximate extinction map. The H-K color excess was adopted as a measure 
of extinction (Lada et al. 1994). This is an effective technique as the spread in H-K colors 
of main-sequence stars is small (~ 0.3 mag). We require 2MASS colors at present, so some 
highly extincted background stars (for which we could in principle estimate Ay from Spitzer 
colors alone) could not be used. But there are enough 2M ASS sources for a map, and the 
accuracy and wide use of this method led us to adopt it. 

We assumed a mean intrinsic color of H-K = 0.175 for the background stars, which 
is the mean value for stars in the Pipe Nebula control field of Lombardi et al. (2006), 
obtained from M. Lombardi (priv. comm.). The conversion to Ay as described above is Ay 
= 16.8E(H-K), where E(H-K) is the color excess. The factor 16.8 differs from the commonly 
used 15.9, due to the adoption of the Indebetouw et al. (2005) near-IR extinction law. Note 
that the maximum uncertainty in intrinsic color translates into an uncertainty in Ay for an 
individual star of 2.5 mag, which is small enough to not affect our results. 

The spatial resolution of the extinction map was set by the need for several stars in 
each spatial bin. This required boxes of size ~ 100", and the map was oversampled by a 
factor 2. The median Ay in each bin was taken to avoid foreground objects and the scatter 
gave the uncertainty. The extinction map is shown in Fig. Id, and interpolated values at 
the positions of the candidate young stars in Table 1. 

The agreement between the total extinction Ay estimates and those towards individual 
young stars is not particularly good. Most of the extinction map Ay values are higher than 
those derived by dereddening individual objects. This is reasonable since many of the objects 
could be near the front of the cloud. Two are lower, which may be explained by localized 
dust. 

The measured peak extinction from the background stars is Ay ~ 45 magnitudes, though 
the true value is likely higher due to the coarse spatial resolution. We estimated the molecular 
hydrogen column density at the peak using N(H2) = 9.4 x 10 20 Ay (Kandori et al. 2005), 
giving N(H 2 ) = 4.2 x 10 22 cm" 2 with a 50 % uncertainty. For the C 18 core in B59, Onishi 
et al. (1999) estimated N(H 2 ) = 1.5 x 10 22 cm -2 from the strength of the C 18 line. 

4.2. Comparison to Gas 



In addition to the CO data of Onishi et al. (1999), data for B59 in CS (J=2-l) and 
N 2 H + (J=l-0) have been obtained at FCRAO with an effective beam size of 60" (C. DeVries, 
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in preparation). CS integrated intensity contours are shown in Fig. le. Source 10 lies close 
to the CS molecular gas peak. The N 2 H + emission (not shown) is centered near this area, but 
appears extended NE-SW, i.e. similar to the extinction. No strong, well-organized outflow 
or infall signatures were detected for either species in these data. 

The positions of the CO outflows within our mapped region from Onishi et al. (1999) 
are marked in Fig. le. These are the positions of the peak contours in the red and blue 
channels defined by those authors. However the large beam size of these observations (2.7') 
makes detailed comparison to the Spitzer sources difficult. It appears that some of the CO 
outflows are aligned with source 11, rather than source 10 which is closer to the extinction 
peak (see Sec. 4.4). 

The total gas mass in the B59 core estimated from the C 18 line is ~ 20 M & within a 
radius of ~ 0.15 pc (Onishi et al. 1999, both adjusted for the different assumed distances). 

The total gas mass in the core can also be estimated from the dust extinction of back- 
ground stars under the assumption that dust within the core provides most of the extinction. 
The conversion from Ay to N(H 2 ) for typical interstellar dust is given above. For a standard 
cloud composition with mean atomic weight per H atom, fi = 1.37 (Lombardi et al. 2006), 
the total core mass integrated within the confines of the C 18 core is M core = 25 ± 1 M , 
where the uncertainty does not include the distance uncertainty. 

These two estimates of core gas mass are in good agreement, considering the different 
spatial resolutions. Some freeze-out of CO is probable in the dense regions. Although solid 
CO has not yet been observed toward B59, sources near the center of the core have deep 
solid H 2 and C0 2 bands (A. Boogert, in prep.). 

4.3. Class 0/1 and I Sources and Bolometric Luminosities 



Based on their SED's there are two deeply embedded objects near the center of the 
B59 core. Source 11 = 2MASSJ171 12317-2724315 (hereafter 2MASS171123) lies ~15" E 
of the quoted position of the 1.3 mm source B59-MMS1 observed with a 22" HPBW beam 
(Reipurth et al., 1996) and the two are probably the same. The total 1.3mm flux integrated 
to the 8o" contour was 725 mJy It was recently detected at 350 /im with the SHARCII 
bolometer array at the Caltech Submilllimeter Observatory (Wu et al. 2006, submitted to 
ApJ). The flux in a 40" diameter aperture is 45.2 ± 8.0 Jy. 

With most of the SED for this source in hand (Fig. 4), we estimated the bolometric 
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luminosity, L bo i= 2.2 ± 0.3 L Q , the bolometric temperature, T bo i= 70 ± 10 K (Myers & 
Ladd 1993), and the ratio of sub- mm (> 350 fim) to bolometric luminosity (L smin / L bo \— 
0.03 ± 0.01). The uncertainty in luminosity does not include the distance uncertainty. 

In addition to the Class I-III spectral classes for young stars discussed above, a fourth 
class, Class 0, is used to define cold sources. Two common criteria for a Class source are 
that T bo i< 70 K and L smm /L bol > 0.5 % (Andre et al. 2000). 

2MASS171123 lies somewhere near the border of this definition, so may be described as 
a Class 0/1 source. 

As noted by Young et al. (2004) and others, heating of dust by the interstellar radiation 
field may contribute to the fluxes of low-luminosity protostars longward of ~ 100 jum. Thus 
high values of L smm /L bo i may be common for low-luminosity protostars at some point during 
their lifetime, and some revision of this criterion may be needed. 

The other nominal Class I object, source 10, is close to the gas and dust extinction 
peaks. It lacks a complete SED at present, but appears spectrally similar to source 11 in the 
mid-IR, so we estimate a bolometric luminosity by scaling, L bo i~ 0.6 L Q , but this is only a 
crude estimate. We will leave its designation as Class I for the present time, pending mm 
and sub-mm observations. 

We are not able to determine directly the final stellar masses for these two Class 0/1 or 
I souces. The models of Myers et al. (1998) suggest that source 11 and source 10 will evolve 
to form stars with 0.5 and 0.3 M , and have ages of t ~ 0.1 — 0.2 Myr. 

It is probable that one or both of the embedded sources near the center of the core we 
detect with Spitzer are the youngest objects and drive the CO outflows seen by Onishi et 
al. (1999). Further study of these objects and their possible disks and outflows with IR and 
mm and sub-mm interferometry is needed. 

Bolometric luminosities (Table 2) for the rest of the candidate young stars were esti- 
mated as follows. For the Class II sources with a < —1.6, late-type photospheres were fitted 
to the dereddened J-band fluxes. The IR excesses add only marginally to the bolometric 
luminosity. For the Class II sources with a > —1.6, we used the absolute J band flux vs. 
bolometric luminosity relation from Greene et al. (1994). For the Flat sources, we have no 
good estimate; only lower limits were obtained by integration under the observed fluxes. 

Estimated bolometric luminosities for the Class II sources range from 0.1-1.8 L Q . The 
central stars provide the bulk of the contribution to the bolometric luminosity for these 
objects. Stellar evolutionary models (e.g. D'Antona & Mazzitelli 1994; Siess et al. 2000) 
show that for plausible ages of t = 0.5 — 1.0 Myr, the stellar masses are ~ 0.1 — 0.9 M , 
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corresponding to mid G to late M dwarfs. There are no clear brown dwarf candidates in the 
present sample. 

4.4. Extended Structure around 2MASS171123 



The Class 0/1 source 2MASS171123 has extended structure at 3.6 and 4.5 /j,m (Fig. 5), 
extending out > 30" (0.02 pc) to the NE and SE. There is a trace of extension at 5.8 /iin 
also. These arcs may trace the edges of an outflow cavity in reflected light. The morphology 
suggests a disk inclined to the line of sight. A possible third extension to the W is confused 
with a superposed background source. 

The apparent outflow axis defined by the E lobes is ~ 50° E of N. This lines up with 
some of the CO outflow peaks (Onishi et al. 1999; see Fig. le). But the low spatial resolution 
of the CO data makes a firm identification of specific outflows with this source difficult. 

The arcs are unlikely to be due to H 2 emission, because the Spitzer IRS spectrum of 
the object shows no evidence of H 2 emission lines at longer wavelengths (A. Boogert, in 
preparation) . 

4.5. Extended Emission 



The Spitzer images of B59 show two kinds of extended emission above the background 
zodiacal emission. 

First is the red (8.0 fj,m) emission seen most clearly in the NW area of Fig. 1/, roughly 
0.8 MJy/sr above background. This is due to galactic cirrus dust emission. The same 
structure is seen in the IRAC 5.8 fim band near the detection limit (not shown). The band 
ratio is I(5.8)/I(8) ~ 0.25. The emission at 24 /iin also broadly traces the same area, with 
1(24) /I(8) ~ 0.4. These flux ratios are in the range expected for galactic cirrus (e.g. Reach et 
al. 2004), with the 8 /zm flux due mostly to the emission in the strong aromatic hydrocarbon 
bands at 7.7 and 8.6 /xm. 

The second extended structure is more interesting: a wide greenish glow centered on 
the core, due to emission in the IRAC 4.5 /im band with a contribution from the blue 3.6 
fim band. This emission is just visible in Fig. If (it may be easier to see in the online 
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version). The band ratio is roughly 1(3.6) /I(4.5) = 0.75. This emission is not detected at 
longer wavelengths. The coincidence of the emission with the cloud core is an indication 
that this is due to scattered galactic light, sometimes called cloudshine (Foster & Goodman 
2006). 

Cloudshine has been seen in many dense cloud cores at visible (Witt & Stephens 1974) 
and near-IR wavelengths (Lehtinen & Mattila 1996; Nakajima et al. 2003; Foster & Goodman 
2006). The JHK data of Foster & Goodman (2006) have been used to study dust structure 
on sub-arcsecond scales (Padoan et al. 2006). It is reasonable to expect cloudshine at longer 
wavelengths as well, though typical interstellar grains will have decreased scattering efficiency 
at longer wavelengths. Longer wavelengths will in principle penetrate deeper into the cloud, 
but the visibility will depend on the sensitivity of the observations and the strength of the 
ambient galactic light. 

This explanation is favored over H2 line emission because of the lack of an obvious jet 
or shock morphology, and because no H 2 lines are seen in the IRS spectrum of source 10 (A. 
Boogert, in preparation). 

5. Discussion 



The Spitzer data reveal the dark cloud core B59 to be a region of active low-mass star 
formation containing a small cluster of young stars at different evolutionary stages. Some 
authors (e.g. Lada & Lada 2003) prefer to reserve the word cluster for groups of > 35 
objects, but we use the term in a general sense here. Some notable aspects and limitations 
of the Spitzer observations are discussed in this section. 

Completeness for Weak Excesses- The sensitivity of this analysis to young stars with 
weak dust excesses is limited by the 24 /xm flux cutoff of 3 mJy A hypothetical young star 
at Ay = 10 with a factor 4 excess at 24 /im would not make the cut if it's unextincted 
photosphere had a flux below 1 mJy. For a main-sequence object at d — 130 pc, this 
corresponds to a mid-K dwarf. For a more likely age of t ~ 1 Myr, this corresponds to an 
object that will become roughly a mid-M dwarf (Baraffe et al. 2002). Later spectral types 
make the cut if they have more substantial excesses. So the sample defined here may be 
incomplete in Class II for late spectral types with weak excesses. 

IRAC Color- Color Plot- Many Spitzer studies of young stars are possible with IRAC 
alone, e.g. Hartmann et al. (2005). Fig. 6 shows an IRAC color-color plot for the candidate 
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young stars in B59 for comparison to IRAC studies of other star-forming regions. The 
colors plotted are as observed, though the IR spectral classes indicated are after correction 
for extinction, where possible. Table 2 has the IR spectral classes before correction for 
extinction. The adopted magnitude zero points were 278 Jy (3.6 /im), 180 Jy (4.5 /iin), 117 
Jy (5.8 /im), and 63.1 Jy (8.0 /im). 

Disk Properties and Age- We have as yet no direct measure of age (e.g. Li absorption 
lines) for any of the candidate young stars. The IR spectral classes give a rough measure of 
age, with Class I being the youngest and Class III the oldest (Kenyon & Hartmann 1995). 

Many authors use a two-class definition for the IR spectral classes: a > for Class I, 
a < for Class II. For comparison in the next paragraph, B59 has a two-class Class I to 
Class II ratio of 3:17, or 18 %. 

The two-class Class I/Class II ratio of 18 % can be compared to Spitzer results on other 
young clusters: 17-34 % for the Tr 37 globule (t ~ 1 Myr, Sicilia-Aguilar et al. (2006); 
Reach et al. (2004)); 25 % for NGC 7129 (t ~ 1 Myr, Muzerolle et al. (2004)); and 33 % for 
the Spokes cluster in NGC 2264 (t ~ 0.5 Myr, Teixeira et al. (2006)). For two clusters in 
Perseus, Jorgensen et al. (2006) find: 29 % for NGC1333 (t < 1 My) and 15 % for the older 
IC348 (t ~ 2 Myr). 

The above comparisons should be taken only as suggestive due to several factors: some 
studies do no extinction corrections; there are different sensitivities for different young star 
selection techniques; and sensitivities differ because of different cloud distances. 

We can also use an argument based on the predicted lifetimes of embedded and T Tauri 
phases (Myers et al. 1987; Kenyon et al. 1990) to estimate the onset of star-formation in 
B59. For a fixed star- formation rate of objects which evolve from Class I {tuf e ~ 10 5 yr) 
to Class II (t/i/ e ~ 10 6 yr) the relative numbers suggest an age of ~ 0.7 x 10 6 yr. Clearly 
the number of objects in each bin in B59 makes this argument only approximate. And we 
have neglected potentially important factors, e.g. the effect of the stellar environment on 
disk lifetimes. 

Multiplicity also complicates the interpretation. The suspected young star pair resolved 
with Spitzer, LkHa 346 NW and SE (itself a probable double), have very different SEDs (our 
IR spectral classes are Flat and Class II, respectively). If coeval, the simple model sketched 
above is wrong. In addition, the frequency of disks may be a function of spectral type, as 
indicated by a recent Spitzer study of IC348 (Lada et al. 2006). 

In Fig. 3, we compared the SEDs of the B59 Class II objects with a > —1.6 to a median 
of a sample of mostly Class II objects in Taurus (D'Alessio et al. 1999). The Taurus SED 
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can be reasonably fitted by optically thick disks with flaring at large radii (D'Alessio et al. 
1999). With time, disks are expected to clear from inside-out due to dust coagulation into 
planetesimals and dynamical clearing. Evidence for such transition disks has recently been 
shown for IC 348 (t ~ 2 Myr, Lada et al.(2006)) and Tr 37 (t ~ 4 Myr, Sicilia-Aguilar et al. 
(2006)). The fact that most of the SEDs of the Class II sources with a > —1.6 in B59 are 
similar to the Taurus SED suggests a comparable age for the cloud, t ~ 0.7 Myr (Kenyon & 
Hartmann 1995). 

We see no clear "transition objects", such as CoKu Tau/4 (Hartmann et al. 2005), those 
which have very little excess in the IRAC bands, then a large 24 /im excess. 

These comparisons taken together suggest an age for B59 of between 0.5 and 1 Myr. 
Deep optical and near-IR photometric and spectroscopic surveys of the cloud are needed to 
ensure the membership and measure more accurate age indicators. 

Star Formation in B59- To date, B59 is the only known star-forming cloud in the Pipe 
Nebula. Onishi et al. (1999) suggest that this is due to its proximity to the B0 star r Oph, 
with a stellar wind having triggered star formation. 

It is notable that the star formation efficiency of B59 is relatively high. Taking just the 
objects near the core, there are 13 young stars in ~ 25 M & of dense gas, or 0.52 M _1 . Of 
the 179 C 18 cores in nearby star-forming regions covered in the NANTEN survey, only a 
few have a comparable efficiency (Tachihara et al. 2002, and references therein). By way 
of contrast, the Taurus cloud as a whole has roughly 200 young stars in ~ 10 4 M of gas 
(Kenyon & Hartmann 1995; Lada et al. 1993). 

Assuming an average mass for the young stars in B59 of 0.5 M , the star formation effi- 
ciency as defined by Lada & Lada (2003), M star / (M star +M gas ) = 0.21, close to the maximum 
of ~ 0.3 seen in their survey of young stellar clusters. 

The surface density of young stars in the inner core region of B59 is ~ 200 pc~ 2 . This 
is a lower limit to the young star population as we require infrared excess for inclusion. The 
surface density is greater than those of 90 % of the young clusters tabulated in Lada & Lada 
(2003), though not as high as much more massive clusters such as S106 and the Trapezium. 
However, the B59 core is smaller than most of the Lada & Lada sample. Comparison to 
the NANTEN core sample cited above again reveals only a few with comparable young-star 
surface density. Recent c2d observations of another small cluster, L1251B, give a comparable 
stellar density to B59 (Lee et al. 2006). 

Compared to all of the cores in the NANTEN CO survey (Tachihara et al. 2002), B59 
is notable for its high gas column density, higher than nearly all of the starless and single- 
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star- forming cores, but similar to their cluster- forming cores (> 10 stars), consistent with 
our survey findings of a small cluster in formation. 

As noted by Lada & Lada (2003), our inventory of young, embedded star clusters, even 
within a few hundred parsecs, is still incomplete. Data from the Spitzer Space Telescope are 
helping to complete that inventory. 

6. Summary 



1) At least 20 candidate young stars in the B59 core have been identified from infrared ex- 
cesses in Spitzer 3.6-70 /zm observations. Spectral energy distributions including corrections 
for extinction are presented. They range from small excesses above photospheric levels to 
rising in the mid-infrared. The SEDs of most of the Class II spectral types with spectral 
index a > —1.6 are similar to the SED of the Class II sample in the Taurus cloud of D'Alessio 
et al. (1999). 

2) Thirteen of the young star candidates are located within ~ 0.1 pc in projection of the 
molecular gas and dust extinction peak. Two Class 0/1 or I sources are identified. One, 
2MASSJ17112317, probably associated with the millimeter source B59-MMS1, with a bolo- 
metric luminosity Lboi ~ 2 L Q , has extended structure at 3.6 and 4.5 /im, possibly tracing 
the edges of an outflow cavity. The other is near the molecular gas and extinction peak and 
previously unrecognized. 

3) The measured extinction through the central part of the B59 core is Ay > 45 mag. There 
are extended infrared emission regions due to galactic cirrus emission and scattered galactic 
light. 

4) Although further study of the young star population in B59 is needed, the data to date 
suggest the B59 core is producing young stars with high efficiency. 
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molecular line data available prior to publication, and R. Gutermuth (CfA) for display 
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Fig 1 - Multiwavelength views of the B59 region, a) 3-color composite in IRAC filters: 3.6 
jum (blue), 4.5 /im (green), and 8.0 jum (red). 6j MIPS 24 /iin image revealing the young star 
cluster. Nonstar-like features seen ±2' roughly N and S of bright objects are artifacts due 
to latent images, c) Source ID key for candidate young stars from Table 1, with IRAC 4.5 
/im image in background. Young star IR spectral classes have the following symbols: Class I 
(circle), Flat (square), and Class II (diamond). Scale bar is for adopted distance of d = 130 
pc. d) Estimated Ay from 2MASS data with Digital Sky Survey Red image. Extinction 
map has resolution ~ 100" and contours at Ay = 10, 14, 19, 30, and 45 mag. e) MIPS 24 
/im image with CS (J=2-l) integrated intensity contours with a 60" beam from C. DeVries 
(priv. comm. ). Contours are at 0.36, 0.50, 0.70, and 0.85 K kms^ 1 . CO red outflow (dark 
x), and CO blue outflow peaks (dark crosses) are from Onishi et al. (1999). These are the 
peak positions of the integrated line flux in red and blue channels. A third blue CO outflow 
peak is off the map to the SE. f) Same as (a) but stretched to show extended emission. 
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0.0 0.5 L0 1.5 2.0 2.5 3.0 

H-K s 

Fig. 2. — 2MASS JHK S fluxes of candidate young stars in B59 (filled circles) and objects 
with J flux > 2.5 mJy and classified as stars by c2d (small x's). Main sequence and giants 
from Bessell & Brett (1988) are the solid lines, and the Taurus T Tauri main locus from 
Meyer et al. (1997) is the dashed line. Both are converted to the 2MASS system using 
Carpenter (2000). The reddening arrow shows the adopted extinction law (Indebetouw et 
al., 2005) converted to Ay (see text) and with Ay = 10 mags. 
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Fig. 3. — Spectral energy distributions (SEDs) of the candidate young stars in B59 as 
observed (filled circles) and after correction for extinction as described in text (open circles). 
Error bars, including absolute calibration uncertainties, are shown, but are smaller than the 
symbols. In each object panel, the source number (Table 1) is in the upper right and the 
IR spectral class after correction for extinction, where possible, is in the upper left. The 
two Class I sources are not corrected for extinction. For the Flat and Class II sources with 
a > —1.6, the median Class II SED of Taurus sources (lower right box, D'Alessio et al. 1999) 
scaled at H (1.65 //m), is plotted for comparison (short dashed lines). Class II sources with 
a < —1.6 have a K7 photosphere scaled at J(1.25 //m) (solid lines). 
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Fig. 4. — Spectral energy distributions of the candidate Class 0/1 or I protostars in B59, 
compared to the infrared source in L1014. For 2MASS171123, the SHARCII 350 data 
point is from Wu et al. (2006, submitted to ApJ). SEST 1.3mm data point of MMS1 is from 
Reipurth et al. (1996). See text for details. 
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Fig. 5. — Close-up of the Class 0/1 source 11 = 2MASS171123, closest to center, through 
the 4 IRAC bands: a) 3.6 /im, b) 4.5 /im, c) 5.8 /im, d) 8.0 jum. Image arcs visible at 3.6 
and 4.5 /im may trace an outflow cavity in scattered light. Scale bar plotted in (b) is for 
assumed distance of d — 130 pc. 
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Fig. 6. — IRAC color-color diagram, as observed, of candidate young stars in B59 with source 
numbers from Table 1. Young star IR spectral classes, after correction for extinction (where 
possible) are indicated with the following symbols: Class I (circle), Flat (square), and Class 
II (diamond). The reddening arrow shows the adopted extinction law (Indebetouw et al., 
2005) for Ay « 30 mags (see text). 



Table I. Candidate Young Stars in B59 a 



No. 2MASS Name Other ID RA(J2000) DEC(J2000) 



K s A V (JHK)9 A v (Ma,p) h 



Comments 



1 J17110392-2722551 LkH« 346 NW 17 11 03.91 -27 22 55.2 

2 J17110411-2722593 LkHa 346 SE 17 11 04.12 -27 22 59.3 



3 J17111182-2726547 



4 J17111445-2726543 



5 J17111626-2720287 



17 11 11.82 -27 26 55.0 



17 11 14.45 -27 26 54.4 



17 11 16.27 -27 20 28.) 



6 J17111631-2725144 KW 2 17 11 16.32 -27 25 14.6 



7 J17111726-2725081 IRAS17081-2721 17 11 17.28 -27 25 08.2 



J17111827-2725491 



17 11 18.13 -27 25 49.3 



9 J17112153-2727417 IRAS17082-2724 17 11 21.50 -27 27 42.3 



10 



11 J17112317-2724315 



12 J17112508-2724425 



17 11 22.10 -27 26 2.0 



17 11 23.18 -27 24 31.5 



17 11 25.08 -27 24 42.7 



10.46 
(0.03) 
9.75 
(0.03) 

14.09 
(0.02) 

11.62 
(0.02) 

11.23 
(0.02) 

10.58 
(0.02) 

13.61 
(0.04) 

>18.4 



12.74 
(0.02) 



8.99 



7.76 



16.52 
(0.11) 



(0.03) (0.02) 
8.79 8.05 
(0.04) (0.03) 



12.62 11.76 
(0.03) (0.02) 



10.42 



9.37 



>18.8 >17.8 



9.63 



(0.02) (0.02) 

9.70 9.07 
(0.02) (0.02) 



8.75 



(0.02) (0.02) 

10.82 8.77 

(0.04) (0.03) 

15.12 11.95 

(0.08) (0.03) 

10.57 8.98 

(0.02) (0.02) 



15.08 
(0.14) 



13.61 11.57 
(0.03) (0.02) 



2.6 
(1.7) 

0.0 
(1.5) 

5.6 
(1.4) 

2.7 
(1.3) 

8.3 
(0.8) 

4.2 
(1.1) 

13.0 
(1.8) 

44 

(8) 

8.7 
(1.7) 



10.0 Triple 1 ' 



IRAS17079-2719 



14.7 
(2.2) 



(0.6) 
9.9 
(0.6) 

19.2 
(1.9) 

22.4 
(2-8) 

7.3 
(0.3) 



14.7 Hq c 
(2.3) 

15.7 
(2.5) 

25.3 
(6.2) 

26.9 
(5.0) 

46.1 
(20.0) 

25.7 Prob. MMSl d 
(5.7) 

27.2 
(5.7) 



Table 1 — Continued 



No. 2MASS Name Other ID RA(J2000) DEC(J2000) 



H 



K s Av(JHK)» Av(Map)' 1 Comments 



13 J17112701-2723485 B59-1 17 11 26.95 -27 23 48.4 



11.88 10.14 
(0.04) (0.04) 



9.08 
(0.03) 



7.5 
(1.5) 



24.4 Triple 6 
(2.4) 



14 J17112729-2725283 



17 11 27.06 -27 25 29.5 



13.17 10.65 
(0.02) (0.02) 



9.11 
(0.02) 



13.7 
(1.3) 



33.2 
(8.6) 



15 J17112942-2725367 



17 11 29.31 -27 25 36.3 



13.28 11.74 
(0.02) (0.02) 



10.70 
(0.02) 



5.1 
(1.3) 



23.6 
(3.7) 



16 J17113036-2726292 



17 11 30.29 -27 26 29.3 



11.91 
(0.03) 



10.00 8.89 
(0.02) (0.02) 



9.3 
(1.4) 



20.7 
(2.7) 



17 J17114099-2718368 IRAS17085-2715 17 11 40.99 -27 18 37.0 



18 J17114182-2725477 B59-2 17 11 41.73 -27 25 50.3 



10.64 8.75 
(0.02) (0.04) 



11.3 
(0.2) 



10.4 
(0.2) 



7.54 
(0.02) 

10.1 
(0.2) 



8.2 
(1.4) 

4.1 

(4.3) 



8.4 
(0.8) 

11.7 Doubled 
(0.6) 



19 J17114315-2730584 



17 11 43.16 -27 30 58.6 



14.09 12.75 
(0.03) (0.03) 



11.91 
(0.03) 



4.1 
(1.4) 



10.1 

(1.2) 



20 J17120020-2720180 LkHa 347 17 12 0.20 -27 20 18.1 



10.54 



9.64 



(0.02) (0.03) 



9.17 
(0.02) 



1.4 
(1.1) 



6.9 
(0.8) 



Note. - 

a Within field covered by both IRAC and MIPS 24 fim. Uncertainties in parentheses. Adopted flux zero points for the 2MASS J, H, and K s filters are 
1594, 1024, and 667 Jy, respectively. 

b LkHa 346 SE is labelled the primary (A) and has a tertiary companion (C) according to Chelli et al. (1995). The AC pair is unresolved by Spitzer. 
LkHa 346 NW, or (B), is the brighter infrared source. 

c Ha emission (Kohoutck & Wchmcycr 2003). 

d Probably MMS1 of Reipurth et al. (1996). See text. 



e The source is triple according to Koresko (2002). The brighter source is the primary (plus tertiary), for which the position and fluxes arc given. 

f The source is double (Rcipurth & Zinnecker 1993), possibly triple (Koresko 2002). The brighter mid-infrared source is the secondary, for which the 
position and fluxes are given. The JHK fluxes are from Koresko (2002). 

g Extinction from dereddening to T Tauri JHK S locus. 

h Cloud extinction from background stars. 



Table 2. Spitzer Fluxes (mJy) of Candidate Young Stars in B59 a 



No. 


Other ID 


3.6 


4.5 


5.8 


8.0 


24 


70 


a" 


IR Class c 




IR Class e 


Lboi(L©)-'' 


1 


LkHa 346 NW 


699 
(9) 


881 
(13) 




1440 

(25) 


6320 
(1260) 


8490 
(40) 


-0.22 


Flat 


-0.35 9 


Flat 


>1.6 


2 


LkHa 346 SE 


218 

(3) 


180 
(2) 


174 

(4) 


141 

(1) 






-1.81 


II 


-1.81 


II 


0.50 
(0.11) 


3 




15.5 
(0.2) 


15.8 
(0.1) 


16.2 
(0.1) 


20.9 
(0.1) 


94.5 
(0.4) 


110 

(5) 


-0.64 


II 


-0.92 


II 


0.14 
(0.08) 


4 




90.1 
(1.1) 


86.6 
(0.7) 


78.4 
(0.5) 


78.1 
(0.4) 


155 
(1) 


210 

(6) 


-1.14 


11 


-1.27 


II 


0.26 
(0.09) 


5 




106 

(1) 


88.0 
(1) 


77.1 
(0.5) 


68.1 
(0.6) 


45.5 
(0.3) 


<100 


-1.64 


II 


-2.04 


II 


0.98 
(0.09) 


6 


KW 2 


211 

(3) 


217 

(3) 


223 
(2) 


228 
(3) 






-0.94 


II 


-1.14 


II 


0.93 
(0.13) 


7 


IRAS17081-2721 


609 

(8) 


711 
(16) 


929 
(11) 


1760 

(24) 


9200 
(1840) 


10100 

(50) 


0.66 


I 


0.02 


Flat 


>2.7 


8 




58.9 
(0.7) 


102 

(2) 


126 
(1) 


169 

(2) 


412 

(3) 


1440 
(10) 


1.11 


I 


-1.04 


11 


1.4 
(1.0) 


9 


IRAS17082-2724 


310 
(4) 


423 
(6) 


601 

(6) 


812 

(9) 


1540 
(17) 


2070 
(10) 


0.21 


Flat 


-0.22 


Flat 


>0.86 


10 




3.2 
(0.1) 


12.5 
(0.1) 


21.6 
(0.2) 


39.9 
(0.2) 


1710 

(29) 


9190 
(30) 


2.16 9 


I 






0.63 
(0.13) 


11 




6.4 
(0.2) 


34.3 
(0.8) 


83.0 
(0.8) 


158.0 
(1.8) 


4715 
(940) 


42100 
(8400) 


3.29 


0/1 






2.2 
(0.3) 


12 




34.2 
(0.3) 


46.4 
(0.4) 


46.2 
(0.3) 


48.1 
(0.3) 


88.7 
(0.5) 




-0.13 


Flat 


-0.85 


II 


0.16 
(0.09) 



Table 2 — Continued 



No. Other ID 3.6 4.5 5.8 8.0 24 70 a b IR Class c a d IR Class e L bo i(L ) / 



13 


B59-1 


146 
(2) 


178 

(2) 




235 
(3) 


831 

(5) 


1780 
(10) 


-0.68 


II 


-1.04 


II 


0.53 
(0.13) 


14 




222 

(3) 


259 
(3) 


250 

(3) 


224 

(3) 


455 
(3) 


710 

(7) 


-0.69 


II 


-1.36 


II 


0.81 
(0.16) 


15 




37.5 
(0.3) 


36.3 
(0.3) 


36.2 
(0.3) 


38.3 
(0.3) 


104 
(1) 


246 
(4) 


-0.93 


11 


-1.18 


II 


0.17 

(0.09) 


16 




243 
(3) 


244 

(3) 


254 
(2) 


253 
(3) 


150 

(1) 


1310 

(6) 


-0.76 


11 


-1.21 


11 


0.74 
(0.16) 


17 


IRAS17085-2715 


471 

(7) 


412 

(6) 


450 
(4) 


670 

(8) 


1050 

(9) 




-0.97 


11 


-1.37 


II 


1.81 
(0.37) 


18 


B59-2 


74.1 
(0.8) 


68.2 
(0.8) 


51.2 
(0.3) 


46.9 
(0.4) 


56.5 
(0.4) 


<100 


-1.57 h 


II 


-1.67 fc 


II 


0.1-0.9 


19 




14.0 
(0.1) 


13.6 
(0.2) 


12.3 
(0.1) 


11.6 
(0.1) 


21.6 
(0.2) 


<100 


-0.99 


11 


-1.19 


II 


0.13 
(0.08) 


20 


LkHo 347 


119 

(2) 


114 

(2) 


101 

(1) 


120 
(1) 


213 
(1) 




-1.13 


II 


-1.20 


II 


0.36 
(0.10) 



Note. — 

a Within field covered by both IRAC and MIPS 24 (im. la errors (in parentheses) arc from photometric fitting and do not include 
absoulute flux error (see text). A blank entry indicates no useful value due to confusion with nearby source (s), or due to lack of data. 

b Spectral index cilog(AFx)/cilog(A) with endpoints at 2.2 and 8.0 £tm uncorrected for extinction, unless different endpoint wavelengths 
in fim are noted. 



c IR spectral class. 



d Same as b after correction for extinction. 
e Same as c after correction for extinction. 

f Bolometric luminosity (or lower limit) estimated as discussed in text. Uncertainties (in parentheses) do not include uncertainty 
distance. 

g On border of Flat and Class II, but longer wavelengths indicate Flat. 
h Spectral index is from 3.6 to 8.0 ^m. 



